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ABSTRACT 

After nearly a decade of quiescence, the soft gamma-ray repeater SGR 1627^1 reactivated on 2008 May 28 
with a bursting episode followed by a slowly decaying enhancement of its persistent emission. To search for the 
still unknown spin period of this SGR taking advantage of its high flux state, we performed on 2008 September 
27-28 a 120 ks long X-ray observation with the XMM-Newton sateflite. Pulsations with P = 2.594578(6) s 
were detected at a higher than 6(t confidence level, with a double-peaked pulse profile. The pulsed fraction in 
the 2-12 keV range is 19% ± 3% and 24% ± 3% for the fundamental and the second harmonic, respectively. 
The observed 2-10 keV flux is 3.4 x 10"'^ erg cm~^ s~', still a factor of 5 above the quiescent pre-burst- 
activation level, and the spectrum is well fitted by an absorbed power law plus blackbody model (photon index 
r ~ 0.6, blackbody temperature kT ~ 0.5 keV, and absorption A^h « 1 .2 x 10^^ cm~^). We also detected a shell 
of diffuse soft X-ray emission which is likely associated with the young supernova remnant G337.0-0.1. 
Subject headings: ISM: individual (G337. 0-0.1) — stars: neutron — supernova remnants — X-rays: individ- 
ual (SGR 1627-41) — X-rays: stars 



1. INTRODUCTION 

Two classes of young isolated neutron stars, the anoma- 
lous X-ray pulsars (AXPs) and the soft gamma-ray repeaters 
(SGRs) are currently believed to host magnetars, i.e. iso- 
lated neutron stars endowed with ultra-strong magnetic fields, 
B - 10''*-10'^ G on the surface and one (or more) order of 
magnitude higher in the interior (Duncan & Thompson 1992; 
Thompson & Duncan 1993, 1995). SGRs and AXPs share 
a number of properties, including rotation periods in the 2- 
12 s range, large period derivatives (P ~ 10~'^-10~^ s s~^), 
inferred surface magnetic dipole field strength B > lO'^ G, 
large and variable X-ray luminosities (exceeding that avail- 
able from the braking of their rotation), and the emission of 
short bursts (see Woods & Thompson 2006 and Mereghetti 
2008 for recent reviews). 

The five known Local Group SGRs were discovered as 
sources of short (<1 s) and intense (lO'^'-lO^^ erg s~') bursts 
of gamma rays that they emit during sporadic periods of ac- 
tivity interrupted by long stretches of quiescence. Once local- 
ized through their bursts and flares, all of them were found to 
be persistent X-ray emitters, with luminosities of about lO-'^- 
10^^ erg s"'. In four SGRs, SGR 1627-41 being the excep- 
tion, pulsations were detected with periods from 5 to 8 s and 
period derivatives of 10"' ^-10"^ s s~'. 

SGR 1627—41 was discovered in 1998 June by the Comp- 
ton Gamma Ray Observatory thanks to the intense bursts it 
emitted at that time, about 100 in six weeks (Woods et al. 
1999). Soon after the discovery of the bursts, its persistent X- 

Based on observations obtained with XMM-Newton, an ESA science 
mission with instruments and contributions directly funded by ESA Mem- 
ber States and NASA. 



ray emission was detected by BeppoSAX at a flux level of '--^ 
7 X 10~'^ erg cm~^ s~' (unabsorbed, 2-10 keV), correspond- 
ing to a luminosity of ^ 10^^ erg s~' for a source distance 
li = 1 1 kpc (Corbel et al. 1999). This led to the possible asso- 
ciation of the SGR with the young supernova remnant (SNR) 
G337.0-0.1 in the CTB 33 complex (Woods et al. 1999; Hur- 
ley et al. 1999). In the following 10 years no further bursting 
activity was reported while the persistent emission steadily 
decayed to about 10^^ erg s~' (for d — 11 kpc; Corbel et al. 
1999), the lowest value ever observed from an SGR. At the 
same time the spectrum softened (Kouveliotou et al. 2003; 
Mereghetti et al. 2006). 

The long-term fading of SGR 1627-41 was suddenly inter- 
rupted by its burst reactivation in 2008 May. This episode was 
associated with a temporary enhancement of the persistent 
X-ray flux and a marked spectral hardening (Esposito et al. 
2008). 

The detection of SGR-Uke bursts from SGR 1627^1 made 
it a bona fide member of the SGR class. However, two strong 
pieces of evidence in favor of this identification, the measure 
of the source period and period derivative, were still missing.^ 
In order to search for pulsations taking advantage of the en- 
hanced flux after the reactivation, we requested a Target of 
Opportunity XMM-Newton observation, that was performed 
as soon as the satellite visibility constraints allowed it. 

2. OBSERVATION, DATA ANALYSIS, AND RESULTS 

Our deep XMM-Newton observation of SGR 1627-41 
started on 2008 September 25 at 00:5 1 :02 UT and lasted about 

' Woods et al. (1999) reported a candidate periodicity of mai'ginal signifi- 
cance at 6.41 s, but the signal was not confirmed by subsequent observations. 
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Fig. 1. — Z| periodogram computed for the pn data of SGR 1627^1. The 
2.594578s pulsed signal and its first harmonic stand out well above the noise 
level. 

120 ks, occupying the whole satelUte revolution number 1611. 

The data were collected with the EPIC instrument, which con- 
sists of two MOS (Turner et al. 2001) and one pn (Striider 
et al. 2001) CCD cameras sensitive to photons with energies 
between 0.1 and 15 keV. The EPIC pn was operated in fuU 
frame mode (time resolution 73 ms) while both MOS units 
were in small window mode (time resolution 0.3 s). All the 
detectors mounted the thick optical blocking filter 

The raw observation data files (ODF) were processed with 
the XMM-Newton Science Analysis Software'^ (SAS version 
8.0.0) and the calibration files released in 2007 August us- 
ing standard pipeline tasks (epproc and emproc for the pn and 
MOS, respectively). The source photons for the timing and 
spectral analysis were accumulated from a circular region (30 
arcsec radius) centered on SGR 1627^1. The background 
events were extracted from source-free regions of the same 
chip as the source. We selected events with pattern 0-4 and 
pattern 0-12 for the pn and the MOS, respectively. For the 
timing analysis, the photon arrival times were converted to the 
solar system barycenter using the SAS task barycen. Photons 
having energies below 2 keV and above 12 keV were ignored, 
owing to the very few counts from SGR 1627-41. A total of 
about 2650 ± 60 counts above the background were collected 
from SGR 1 627-4 1 by the pn detector in the 2-12 keV range, 
760 ± 40 by the MOS 1 detector, and 900 ± 40 by the MOS 2 
detector. 

2.1. Timing analysis 

For the timing analysis we started with the high time- 
resolution pn data. We searched for pulsations in the 0.15- 
15 s range^' using the test (Buccheri et al. 1983), with the 
number of harmonics n being varied from 1 to 5. The most 
significant peak occurred at about 2.6 s for n = 2 (see Figure 
1). Taking into account the number of searched harmonics 
and periods (1,647,257), the value of 80.72 for this peak 
corresponds to a chance probability of ^ 9 x 10"'^ (that is a 
>6ct detection). 

We refined the period measurement by folding the pn plus 
MOS data at the best period inferred through the test and 
studying the phase evolution by means of a phase- fitting tech- 

See http://xmm.esac.esa.int/extemal/xmm_data_analysis/ 
' ' The frame time of the pn in full frame mode implies a Nyquist limit of 
about 146.6 ms. 
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Fig. 2.— Epoch-folded pulse profile of SGR 1627^1 (2-12 keV) for the 
pn (top) and MOS data (bottom). The arrival time of the pulse minima (arbi- 
trarely adopted as phase 0) was 54734.0000000(4) MJD (TDB). 

nique (see Dall'Osso et al. 2003 for more details). The re- 
sulting best-fit period was P = 2.594578(6) s (epoch 54734.0 
MJD; hereafter all uncertainties are at la confidence level, 
unless otherwise specified). The current data set a 3 u upper 
hmit on the period derivative of |P| < 6 x 10"'" s s~' . 

The background-subtracted lightcurves folded at the best 
period are shown in Figure 2. The pulse profile is double 
peaked, with a pulsed fraction (semiamphtude of sinusoidal 
modulation divided by the mean source count rate; see Israel 
& Stella 1996 for more details) of 19% ± 3% and 24% ± 3% 
for the fundamental and the second harmonic, respectively 
(the above values correspond to 16% ±4% in terms of root 
mean square, consistently with the 10% upper hmit reported 
in Woods & Thompson 2006). Within the statistical uncer- 
tainties the pulsed fraction is energy independent and no sig- 
nificant pulse shape variations as a function of energy were 
found by dividing the counts in soft and hard energy inter- 
vals. To assess this, we compared the folded lightcurves in 
various energy ranges using a two-dimensional Kolmogorov- 
Smirnov test (Peacock 1983; Fasano & Franceschini 1987). 
The results show that the probability that they come from the 
same underlying distribution is always larger than 30%. 

Knowing the spin period, we went back to archival data and 
looked for pulsations with enhanced sensitivity. However, our 
searches in the previous XMM-Newton observations, as well 
as in BeppoSAX, ASCA and RossiXTE data,'^ were inconclu- 
sive. The most constraining limits on the pulsed fraction were 
obtained from BeppoSAX (1998 August) and ASCA (1999 
February) archival observations of SGR 1627^1. A Fourier 
analysis of the background-subtracted light curves was per- 
formed by means of the method described in Israel & Stella 

Public available Chandra data of SGR 1627-41 could not be used owing 
to their inadequate frame times. 
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TABLE 1 

Spectral results for SGR 1 627-4 1 . 



Parameter 




Value 






BB 


PL 


PL+BB 


Nn (1022 cm-2) 

r 

kT (keV) 
R(km) 


1.7 ±0.4 

2 'i-^-2 
0.041 ±0.003 


5 2"*^ * 

"•'^-0.08 


12±2 
1.2±0.1 


X?/dof 


1.38/152 


1.19/152 


1.07/150 


Flux" (10"" erg cm-2 ^-l^ 


3.3 


3.3 


3.4 



Note. — The abundances used are those of Anders & Grevesse (1989) 
and photoelectric absorption cross-sections from Balucinska-Church & Mc- 
Cammon (1992). Blackbody temperature and radius are calculated at infinity 
and assuming a distance to the source of 1 1 kpc (Corbel et al. 1999). 



"Observed flux in the 2-10 keV range. 

(1996) and no significant periodicity was found in the period 
interval 2.4-2.8 s (identified by considering the 3(t upper limit 
on P inferred above). The 3a upper limits on the pulsed frac- 
tion (as defined above), computed according to Vaughan et al. 
(1994), are about 63% and 51% for the BeppoSAX and A5CA 
data sets, respectively.^^ 

2.2. Spectral analysis 

For the spectral fitting (with XSPEC version 12.4; Arnaud 
1996), data were grouped so as to have at least 30 counts per 
energy bin. The ancillary response files and the spectral redis- 
tribution matrices were generated with the SAS tasks arfgen 
and rmfgen, respectively. We jointly fitted the spectra by pn, 
MOS 1, and MOS 2 to blackbody, power law, and blackbody 
plus power law models, all corrected for interstellar absorp- 
tion (see Table 1 for the best-fit model parameters). 

The data disfavor the blackbody model (x? = 1-38 for 152 
degrees of freedom (dof)) and also the power-law fit yields a 
relatively high value (x^ — 1.19, 152 dof) with structured 
residuals (see Figure 3). A better fit is obtained by a power 
law plus blackbody model (x^ = 1.07, 150 dof). The best-fit 
parameters are photon index T ~ 0.6, blackbody temperature 
kT ~ 0.5 keV, and absorption A^h ~ 1-2 x 10^ cm . For a 
distance of 11 kpc (Corbel et al. 1999) this corresponds to a 
luminosity of ^ 10^'^ erg s"' (2-10 keV, unabsorbed). 

We compared the current spectral shape of SGR 1621 -A\ 
with those inferred from the latest observations by means of 
simultaneous fits. Apart from an overall normalization factor 
(to account for the different luminosities), the current spec- 
trum is consistent with that observed by 5w;/f/XRT during the 
onset of the 2008 May outburst, while it is significantly harder 
than that observed by XMM-Newton before the outburst, in 
2008 February (Esposito et al. 2008). 

2.3. Imaging analysis 

The energy-coded image of the field around SGR 1627^1 
(Figure 4) shows that the source is embedded in a complex 
region of diffuse emission. A shell-shaped structure (^ 2 ar- 
cmin diameter) with a softer spectrum than that of the SGR, as 
well as a bright spot of remarkably hard diffuse emission (^ 2 
arcmin to the southwest of SGR 1627-41) are apparent in Fig- 
ure 4, where red and blue correspond to soft and hard X-ray 
photons, respectively. The combination of large absorption, 

" These limits do not take into account the contribute from the diffuse 
structures described in Section 2.3, that are not resolved in the BeppoSAX 
and A5CA data. 



source confusion and small photon statistics makes spectral 
analysis of such diffuse structures rather difficult. Selecting 
photons from a ~ 100" x 50" rectangular region located be- 
tween SGR 1627-41 and the hard spot, we observe the emis- 
sion of the very faint diffuse shell to peak between 2 and 3 
keV. Using Chandra imaging data of the field, we estimated 
the contribution of unresolved point-like sources to the dif- 
fuse shell to be of order 30% in the 1.7-3.1 keV energy range. 
The brightest point sources in the region are also resolved in 
the EPIC data below 1.7 keV, where the diffuse emission is 
not seen. The spectrum of the shell is adequately described 
(Xr = 1 • 13, 106 dof) by an absorbed (A^h 8 x 10^2 cm'^) hot 
plasma model (pshock in XSPEC, kT keV), with clear in- 
dications for emission lines from S (at 2.46 and 3.10 keV) and 
possibly from Fe (at 6.4 keV) and is therefore consistent with 
emission from a very absorbed SNR. The net (point source 
subtracted) surface brightness of the structure is of ~ 10~^^ 
erg cm~^ arcmin-^ in the 1.7-3.1 keV band. 

The spectrum of the second, harder, diffuse structure was 
extracted using a 20 arcsec radius circle around the apparent 
centroid of the feature. The source emission emerges from 
background above ^ 2 keV and peaks around 5 keV. Signifi- 
cant signal is detected up to 12 keV. A thermal plasma model 
(mekal in XSPEC, abundance parameter left free) yields an 
acceptable description of the data (Xr = 1-27, 53 dof), al- 
though the absorbing column and the temperature are not well 
constrained (A^h > 1 .3 x 10^^ cm"^, kT > 13 keV at 90% con- 
fidence level for a single parameter). The possible presence 
of an emission feature at ^ 6.2 keV makes the fit improve 
after freeing the redshift parameter (x^ = 1 -08, 52 dof; best- 
fit z ~ 0.1, pointing to Fe emission at 6.9 keV in the source 
frame). The 2-10 keV flux of the source within the aperture 
is of 1.2 X 10~'^ erg cm-2 ^-i ^ possible interpretation 
is that this source is a background cluster of galaxies. We 
note that a simple absorbed (A^h ~ 1.4 x 10^^ cm'^) power 
law (F ~ 1 .6) also yields an acceptable fit (x? = 1 -27, 53 dof). 
Thus, the source could also be a pulsar wind nebula, unrelated 
to SGR 1627-41. 

3. DISCUSSION AND CONCLUSIONS 

We have discovered pulsations with period of 2.6 s in 
SGR 1627^1, which was the only known magnetar with 
no spin period yet measured. The pulse shape is double 
peaked, with pulsed fractions of ~ 19% and 24% in the two 
harmonics. No dependence of the pulse profile on the energy 
band has been found. Among the known magnetars, only 
IE 1547.0-5408 has a shorter period (2.07 s; Camilo et al. 
2007) than that of SGR 1627-41. We searched for a similar 
periodicity in all available archival data with sufficient time 
resolution, but no pulsation was found. Relying only on 
the last XMM-Newton observation we could set an upper 
limit on the period derivative, |P| < 6 x 10~^*^ s s~^ which 
is not particularly constraining for a magnetar, and yields 
B < 1.2 X lO'^ G. Based on the peak luminosities of bursts 
observed from SGR 1627-41, its dipole magnetic field was 
estimated to be 5 > 1 .8 x lO^"* G (Woods et al. 1999; Esposito 
et al. 2008). This limit can be used to infer a period derivative 
/" > 1.2 X 10"'^ s s~\ a spin-down luminosity E >8 x 10^^ 
ergs~\ and a characteristic age Tc = P/{2P) < 3.4 kyr. 
Overall, and lacking any information on the pulsation period 
and profile in the past, we can only conclude that the timing 
properties of SGR 1627-41 reported here appear quite typical 
of a member of the SGR class. 
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Fig. 3. — EPIC/pn spectrum of SGR 1627—41. Top: data and best-fit power law plus blackbody model (solid line); the dotted (blackbody) and dashed line 
(power law) show the individual components. Middle: residuals from the power-law best-fit model in units of standard deviation. Bottom: residuals from the 
power law plus blackbody best-fit model in units of standard deviation. The data have been rebinned in the plot to better visualize the trend in the spectral 
residuals. 





Fig. 4. — EPIC/pn image of the field of SGR 1627^1 (1 1' x 7'). Data from the long 2008 September observation are combined with previous XMM-Newton 
data, totaling to 147.8 ks exposure time. Photon energy is color-coded: red coiTesponds to the 1.7-3.1 keV energy range, green to 3.1-5.0 keV, blue to 5.0-8.0 
keV. North is up, east is left. SGR 1627^1, in white, is the brightest source in the field. The image unveils diffuse X-ray emission from SNR G337. 0-0.1, 
prominent in the softest energy range. A patch of hard diffuse emission is also apparent, southwest of the SNR, possibly due to a background, very absorbed 
cluster of galaxies. 



Our observation caught SGR 1627-41 at a flux of 
3.4 X 10"'^ erg cm~^ s~', the highest seen with XMM-Newton 
and a factor of 5 above the level preceding its 2008 May 
burst activation (see Figure 1 of Esposito et al. 2008 for a 
complete flux history of the source). Its spectrum was well 
fitted in earlier data by a single-component model, either 
a heavily absorbed power law or a blackbody (Mereghetti 
et al. 2006). The higher-quality spectrum of the new data 
cannot be adequately fitted by these simple models, while it 
is consistent with a power law plus blackbody model, and it is 
significantly harder than in any observation of SGR 1627-41 



taken before 2008 May (Kouveliotou et al. 2003; Mereghetti 
et al. 2006). The presence of a soft thermal component 
and the spectral hardening as the source moves from a 
quiescent to a (burst) active phase makes SGR 1627-41 akin 
to other SGRs and AXPs (e.g. Mereghetti et al. 2005; Rea 
et al. 2005). This latter behavior can be interpreted in the 
framework of the "twisted-magnetosphere model" as due 
to a progressive growth of the shear in the magnetosphere 
as magnetic helicity is transferred from the internal to the 
external field (Thompson et al. 2002; Fernandez & Thompson 
2007; Nobili et al. 2008). 
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SGR 1627^1 has been proposed to be associated (Woods 
et al. 1999; Hurley et al. 1999, but see also Gaensler et al. 
2001) to the radio SNR G337.0-0.1 (Frail et al. 1996; 
Whiteoak & Green 1996). The shell of diffuse soft radiation 
(Section 2.3), consistent with emission from a young SNR 
at the distance of SGR 1627^1, is Ukely to be the first 
detection of this SNR in X-rays. Due to the complexity of 
the region, the extension of the SNR at radio frequencies is 
debated: Sarma et al. (1997) claimed a diameter of ~ 95 
arc sec with SGR 1627^1 outside the shell, while Brogan 
et al. (2000) proposed a different morphology extending the 
SNR size towards the southwest and encompassing the SGR. 
The position and extension of the diffuse X-ray emission are 
consistent with the latter hypothesis, supporting the possible 
association between SGR 1627^1 and G337.0-0.1. 
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